Stability of Sarma phases in density imbalanced electron-hole bilayer systems 
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We study excitonic condensation in an electron-hole bilayer system with unequal layer densities 
at zero temperature. Using mean-field theory we solve the BCS gap equations numerically and 
investigate the elTects of intra-layer interactions. We analyze the stability of the Sarma phase with 
k, — k pairing by calculating the superfiuid mass density and also by checking the compressibility 
matrix. We find that with bare Coulomb interactions the superfiuid density is always positive in the 
Sarma phase, due to a peculiar momentum structure of the gap function originating from the singular 
behavior of the Coulomb potential at zero momentum and the presence of a sharp Fermi surface. 
Introducing a simple model for screening, we find that the superfiuid density becomes negative in 
some regions of the phase diagram, corresponding to an instability towards a Fulde-Ferrel-Larkin- 
Ovchinnikov (FFLO) type superfiuid phase. Thus, intra-layer interaction and screening together 
can lead to a rich phase diagram in the BCS-BEC crossover regime in electron-hole bilayer systems. 

PACS numbers: 73.21.Ac,03.75.Hh,03.75.Ss 



I. INTRODUCTION 

Recent advances in the trapping and cooling down 
to degeneracy of ultracold Fermi gases have revived in- 
terest in the ground state phases of these systemsii"— 
In a two-component Fermi system with equal densities 
attractive interactions between different species lead to 
Bardeen-Cooper-Schrieffer (BCS) pairing in the weak 
coupling limit and Bose-Einstein condensation (BEC) in 
the strong coupling limit. When the densities are imbal- 
anced more exotic phases are expected to follow, such as 
the Sarma phase^^ with zero center-of-mass momentum 
and Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) phasei^ii^ 
with finite center-of-mass momentum. There is a growing 
literature on the possible phases of two-component Fermi 
gases with population and mass imbalanceJ^^— The ex- 
perimental efforts in ultracold Fermi gases are in their be- 
ginning stage and so far only phase separation between a 
superfiuid and a normal phase has been observed i^iSi^!^ 

Semiconducting electron-hole bilayer systems offer an- 
other realization of a two-component Fermi system 
with which the exotic phases can be studied. Forma- 
tion of excitons between spatially separated electrons 
and holes and their subsequent condensation have long 
been predicted^"^^ and arguably observed nearly 30 
years later experimentally.*^ The phase diagram of sym- 
metric electron-hole bilayer systems (equal mass and 
layer density) is most reliably calculated by quantum 
Monte Carlo simulations Recent success in fabricat- 
ing closely spaced semiconducting electron-hole bilayer 
structures--'— and the ability to control the densities of 
individual layers make the investigation of Sarma and 
FFLO phases very timely. In fact, experiments support- 
ing evidence of a transition from the Fermi liquid phase 
to an excitonic condensate have been recently reported 
through Coulomb drag measurement o^^'^^ in such struc- 



tures. 

The BCS-BEC crossover in an electron-hole bilayer 
system with unequal electron and hole densities was re- 
cently studied in Ref. [H within a BCS mean field ap- 
proach. Sarma and FFLO phases were found to be stable 
in some range of densities; electron-hole bilayers appear 
thus promising candidates for the detection of such elu- 
sive phases. 

In this paper we extend the work of Ref. 5^ by in- 
cluding the in-plane Coulomb interactions, that were ne- 
glected there, as well as some screening effects. We find 
that the effect of intra-layer Fock energy quantitatively 
changes the phase diagram moving the normal-condensed 
phase boundary to lower densities. Comparing energy 
of the condensed phase with that of the normal phase, 
we map out the phase diagram in the average density- 
population polarization plane. We check the "local" 
stability of the Sarma phase with respect to competing 
FFLO order by calculating the superfiuid mass density 
and identify a negative superfiuid mass density with an 
instability towards an FFLO phase. We calculate also 
the compressibility matrix in order to investigate possi- 
ble instabilities towards phase separation. 

Finally, we consider the effect of gate layer screening, 
which proves especially important in the discussion of 
the local stability of the Sarma phase. At zero tempera- 
ture, the simultaneous presence of the singularity in the 
Coulomb potential and of a sharp Fermi surface produces 
in fact a logarithmic divergence in the momentum depen- 
dence of the BCS gap function which makes the Sarma 
phase always locally stable against the FFLO phase. The 
inclusion of some form of screening removes this pecu- 
liar behavior, thus recovering the instability towards the 
FFLO phase in some region of the phase diagram. The 
intra-layer interactions and screening effects give there- 
fore rise to a rich phase diagram in the crossover region 
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between the BCS-like high density state and the BEC of 
low density excitons, showing the possibihty to observe 
exotic superfluid phases as the population polarization is 
changed. 

The rest of this paper is organized as follows. In the 
next section we outline the mean-field theory for electron- 
hole bilayers and provide the set of self-consistent equa- 
tions for the quasiparticle energies and gap function. In 
Section III, after a brief remark about our computational 
procedure, we present our results for the quasiparticle 



properties and phase diagram of the system. We con- 
clude in Section IV with a summary and outlook. 



II. MEAN-FIELD THEORY 

The Hamiltonian describing electrons and holes in a 
bilayer system interacting with the Coulomb potential 
can be written as 
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The basis states for electrons and holes are chosen to 
be plane wave states labeled by two-dimensional wave 
vectors k as is conventional for a uniform system. The 
operators at/aj^ (^k/^k) are creation/annihilation opera- 
tors for electrons (holes) respectively. The single particle 
energies are denoted by e£, and the matrix element [/q 
with respect to plane wave states becomes the Fourier 
transform of the corresponding two-body Coulomb inter- 
action U{y) 
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where U°'°',U^'', and U°'^ denote the electron-electron, 
hole-hole and electron-hole Coulomb interactions, respec- 
tively, A is the area of a layer and d is the inter-layer 
separation. We disregard the spin degrees of freedom. 

The bilayer system is characterized by the electron 
and hole densities, or equivalently by the average den- 
sity parameter (average distance between particles in 
the plane in units of Bohr radius as), the population 
polarization a (characterizing population imbalance in 
terms of the ratio of density difference and total density) 
defined by 
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and the inter-layer separation d. 

The solution of the mean-field Hamiltonian at zero 
temperature (T = 0) is given by the following coupled 



integral equations 
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where ek = (Ck + (^ith ^ h^h^ /2m„ i = a,b), 

the mean chemical potential = (/if, -I- /i6)/2, while 



ASk = A& + ^E^k"k'(/k^-/i) 
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At finite temperature, the occupation functions f^{E^) 
go from the step function to the Fermi-Dirac distribution. 

Given the electron and hole chemical potentials and 
fib, these equations can be solved numerically to obtain 
the unknown functions Ak,^k and Ai?k. Physically, Ak 
is the BCS (s-wave) gap function while E^ are the quasi- 
particle excitation energies in the superfiuid phase. In the 
absence of intra-layer interaction ^k is just the average of 
the free electron and hole dispersions (with respect to the 
corresponding chemical potentials). Intra-layer interac- 
tion modifies the free dispersions by the inclusion of the 
exchange (Fock) interaction, as explicitly considered in 
Eq. ®. 
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For fixed number of particles the chemical potential 
values are adjusted to satisfy the number equations 
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In the above mean-field description of the electron-hole 
bilayer we have used the bare Coulomb interaction given 
in Eq. (2). In realistic systems, the interactions enter- 
ing the model hamiltonian of Eq. (1) should be modified 
to include many-body effects such as exchange and cor- 
relation and external potentials. These effects are de- 
scribed by a screening function which usually decreases 
the strength of the bare Coulomb interaction for elec- 
trons and holes in the normal phase. However, the 2D 
screening due to intra- and intcr-layer interactions is dif- 
ficult to take into account properly for the condensed 
phase<^ In order to see the qualitative effects of screen- 
ing we consider the mechanism of gate screening which 
can be taken into account in a simple way. In this mecha- 
nism the Coulomb potential of a point charge is replaced 
by that of a dipole consisting of the point charge and its 
image behind the metallic gate. We have approximately 
modelled the screening by the gate potential by taking 
the intra- and inter-layer interactions to be 
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respectively, where the parameter k is a screening wave 
number. In recent experiments with metallic gates to 
control the charge densities, the separation between the 
gate and 2D layer is about ^ 250n m^°i^^ . Thus, the im- 
age charge is 500 nm away from the real charge and we 
may assume that for distances larger than 500 nm, the 
Coulomb potential will be screened. In our calculations 
we take the screening length associated with gate screen- 
ing to be ^ 20aB, i.e., k = l/20aB. We have also checked 
other values of k and found that the results are largely 
insensitive in the range 40 > l/(KaB) > 5. 



III. RESULTS AND DISCUSSION 

A. Numerical Procedure 

We solve the gap equations by representing the un- 
known functions on a grid of /c-points (after angular in- 
tegration) and using a non-linear root finding scheme for 
the function values on grid points. For balanced pop- 
ulations an iterative scheme provides a robust method 
of solution. For imbalanced populations we employ a 
root finding scheme for the function values on grid points 



and chemical potential values. We start with the equal 
density solution at the same average density and create 
imbalance first at a small finite temperature and then 
decrease the temperature (using the solution from the 
previous step as input) until results do not change with 
temperature any more. The integrals are evaluated us- 
ing Gaussian quadrature. The finite temperature is nec- 
essary to obtain smooth functions and gradients for the 
Newton-Raphson root finding algorithm. We found it 
necessary to introduce up to three different grids for inte- 
gration to handle "discontinuities" at low temperatures, 
when one type of occupation number becomes equal to 
unity = 1 (A = -I- or — ) in a region of fc-space (and 
causing the integrands to vanish there). 

The so-called Sarma states obtained in this way are of 
the following BCS form 
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where the resulting wave function has a certain range 
of k states (the set denoted by TV) occupied with quasi- 
particles of the BCS theory giving rise to population im- 
balance. The region 7?. is where the quasi-particle energy 
E'^ becomes negative, i.e. less than that of the ground 
pair energy and the corresponding quasi-particle occupa- 
tion becomes unity. Incidentally, the quasi-particles of 
BCS theory are just electron or hole states at that wave 
vector k. Outside the set TZ we have pairs k, — k of elec- 
trons and holes. Therefore, at T = there can be one 
or two Fermi surfaces depending where the set of k e 7?. 
vectors are. These topologically different phases will be 
called Sarma-1 (SI) and Sarma-2 (S2). These states have 
also been called breached pair states^^ displaying a Fermi 
surface together with a condensate. The gap function is 
non-zero but there are gapless excitations. 



B. Quasi-Particle Properties 

In the following we present physical quantities in Ry- 
dberg units, i.e. length is measured in effective (exci- 
tonic) Bohr radius qb = momentum in l/a^ and 

energy in effective Rydberg (Ryd= = 2^)- The 

reduced mass m is defined by 1 /m = 1 / nia + 1 / rrib where 
ma = rue and nib = nih are the band mass of the elec- 
tron and hole, respectively. In the numerical calculations 
we specialize to GaAs system parameters with mass ratio 
ma/rrib — 0.07/0.30 and background dielectric constant 
e = 12.9. 

Representative solutions with one and two Fermi sur- 
faces at r = are illustrated in Fig. [T] (bare Coulomb in- 
teraction) and Fig. [5] (screened interactions) for d = as . 
The figures show the gap function Ak, the quasi-particle 
energies E^ and their average i?k on the left panels, and 
the electron and hole occupation numbers na{k),nb{k) 
on the right panels. At T = in the ground state. 
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the quasi-particle levels with negative energy are occu- 
pied, positive energy levels are empty. The two different 
type of excitation branches are split both due to different 
electron-hole mass and chemical potential values. When 
one of the spectra crosses the zero energy axis, popula- 
tion imbalance is created. If the negative energy region 
includes the origin at fc = 0, the ground state has one 
Fermi surface, otherwise it has two. The two cases are 
denoted by SI and S2 respectively. The top panel in each 
figure shows an S2 phase and the bottom panel shows an 
SI phase. Since the quasi-particle energy branch is con- 
tinuous the system still has gapless excitations. A close 
investigation of the gap function Ak in the absence of 
screening (Fig. [1]) shows that it has a cusp at the zero 
crossings of the quasi-particle energy, corresponding to 
a divergence in the derivative of Ak. This divergence 
has important consequences on the stability of the Sarma 
phase at T = 0, as discussed below. 



C. Superfluid mass density, compressibility matrix 
and the stability of the Sarma phase 

The "local" stability of the Sarma phase with re- 
spect to phases of the FFLO type is usually as- 
sessed by calculating the superfluid mass density (phase 
stiffness).— 1^>^ This quantity should be positive in a 
stable state and a negative value is identified with an 
instability towards an FFLO phaseJ^ii^i^ Clearly, the 
positivity of the superfluid mass density guarantees only 
that the Sarma phase is a local minumum of the energy 
with respect to fluctuations of the gap parameter associ- 
ated with pairing of the FFLO type, and does not exclude 
the possibility that an FFLO phase with finite pair mo- 
mentum can be a global mimimum of the energy. When 
this happens, the local stability of the Sarma phase ac- 
tually corresponds to metastability. 

The superfluid mass density is given by^i 



the bare Coulomb interaction one has 
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where 13 is the inverse temperature. At T 
pression can be written aa^ 
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where k^ are the roots of E^ with A = ±. At zero tem- 
perature the last expression involves the derivative of Afe 
(through the derivative of E^) at the zero crossings of 
E^ . As mentioned above, our calculations are carried 
out at nonzero but small temperature. We have found 
that this derivative diverges logarithmically as T 0. In 
particular, we have demonstrateed analytically that for 
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where k* is the zero crossing point at T = as ^ fc* 
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which we have also checked numerically (Fig. [3]). This di- 
vergence is due to the simultaneous presence of the long 
range Coulomb interaction, which is singular at g = 0, 
and the discontinuity of the Fermi function at T = 0. 
Finite temperature and/or screening effects, smear out 
these singularities thus removing the divergence. The 
presence of this divergence at T = was overlooked in the 
previous mean field study of the imbalanced electron- hole 
bilayer system^. As a matter of fact, making this diver- 
gence emerge from the numerical calculation requires the 
achievement of very low temperatures in the calculations 
and an extreme precision in the numerical integration. 

The presence of this divergence is particularly mean- 
ingful for the analysis of the local stability of the Sarma 
phase at strictly T — 0. The diverging derivative makes 
in fact the negative contribution to ps vanish, thus im- 
plying that for the unscreened Coulomb interaction the 
Sarma phase is always locally stable at T = 0. This re- 
sult is interesting as a matter of principle, as it offers 
an "extreme" example, where the argument by Forbes 
et al.^^ that mass ratio and momentum structure of the 
interactions should favor the stability of Sarma phase is 
completely effective. On the other hand we expect it 
to have few practical consequences, since the stability of 
the Sarma phase induced by this divergence is very frag- 
ile with respect to finite temperature and/or screening 
effects. As Fig. [3] clearly shows, quite small tempera- 
tures suffice to smear out the divergence in the deriva- 
tive. Alternatively, the simple screened interaction makes 
the divergence to disapper even at T = 0, as also shown 
in Fig. [3] As a result, in the presence of screening, the 
Sarma phase will indeed be locally unstable in certain 
regions of the — a plane, as discussed in the next sub- 
section. 

The mechanical stability of the system with respect 
to phase separation requires the compressibility matrix 
dpi/dNj to be positive definite. We have therefore cal- 
culated the compressibility matrix across our phase dia- 
gram to check also this stability. When the intra-layer 
Coulomb interaction is neglected, the compressibility ma- 
trix develops negative eigenvalues across most of our 
phase diagram (restricting the stable region only to small 
values of r^) in agreement with the findings of the re- 
cent work by Yamashita et alM: However, as it was 
already argued in Ref. [H, this apparent dominant insta- 
bility towards phase separation is an artifact occurring 
when the intra-layer Coulomb repulsion is artificially ex- 
cluded from the calculation. It should therefore not be 
taken seriously. In particular, we have verified explic- 
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FIG. 1: (color online) Gap function and quasi-particle energies with bare Coulomb interactions for m,^/m,h ~ 0.07/0.30 and 
d = ub- The upper panel shows a Sarma-2 phase at rs = 3 and a = —0.3 with excess holes (heavy majority species). The 
lower panel shows a Sarma-1 phase at = 5 and a — 0.5 with excess electrons (light majority species). Occupation numbers 
are shown on the right. 



itly that in our calculations with Coulomb intra-layer re- 
pulsion, the Hartree term, which increases linearly with 
the distance da between the metallic gates and the elec- 
tron/hole layers, washes out completely phase separation 
from our phase diagram of Fig. 2] already for distances 
do of the order of 5-10 as, well below the typical gate- 
to-layer distances in current devices. We thus conclude 
that, contrary to what happens in cold atom systems, 
phase separation is not an issue in electron-hole bilayer 
systems. 



D. Phase Diagram at d = as 

In this section we present the phase diagram result- 
ing from the comparison of the energies of the Sarma 
and normal phases and from the stability analysis dis- 
cussed in the previous section. We set the inter-layer 
separation equal to one effective (excitonic) Bohr radius 
d = aB- To make contact with previous literature, we 
present in Fig. |4]the phase diagram for progressively re- 
fined approximations corresponding to the inclusion in 
the calculations of: (i) bare inter-layer interactions only, 
(ii) bare inter- and intra-layer interactions, (iii) screened 
inter-layer interactions only and (iv) screened inter- and 
intra-layer interactions. 

For bare interactions, the superfluid density is always 
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FIG. 2: (color online) Gap function and quasi-particle energies with screened Coulomb interactions for rrie/rah ~ 0.07/0.30 
and d = as- The upper panel shows a Sarma-2 phase at ra — 2.5 and a = 0.2 with excess electrons. The lower panel show a 
Sarma-1 phase at Ts = 5 and a = 0.5 with excess electrons. Occupation numbers are shown on the right. The gap function has 
less variation and the divergence in the derivative at the zero crossings disappears. 



positive and the Sarma phase is "locally" stable, due to 
the mechanism explained in Sec. IIII CI Therefore, in the 
top panel of Fig. |4]we do not show any FFLO phase, but 
our calculations do not rule out the possibility of a first 
order transition to an FFLO phase with a finite FFLO 
modulation momentum q as found in Ref. [STi Two topo- 
logically distinct Sarma phases, Sarma-1 with one Fermi 
surface and Sarma-2 with two Fermi surfaces, are present 
in the phase diagrams. The effect of the intra-layer re- 
pulsive interactions is to favor the normal phase with 
respect to the condensed phases, thus shifting to higher 
vaues of the boundary between normal and condensed 
phases (right panels). The two bottom panels of Fig. |4] 
presents the phase diagram when the gate screening is 
taken into account. With inter-layer interactions only, 
the Sarma phase becomes unstable for a large portion 



of the phase diagram, especially with excess holes, i.e. 
a < 0. (bottom left panel) There is no stable S2 phase 
in this case. Switching on the intra-layer interactions re- 
duces the space occupied by the FFLO phase in our phase 
diagram and restores the S2 phase in some region of the 
phase diagram. This result is physically quite sensible, as 
the FFLO modulations of the gap parameter should be 
unavoidably accompanied by some modulations of the 
density in real space. In the presence of the Coulomb 
intra-layer repulsion such density modulations are ener- 
getically expensive, thus hindering the FFLO phase with 
respect to the Sarma phasei^ 

A quite rich phase diagram is therefore obtained when 
both intra-layer and screening effects are present. The 
presence of locally stable Sarma phases confirms the ex- 
pectation that isotropic translationally invariant gapless 
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FIG. 3: The derivative of the gap function dAk/dfc as a function of k at various temperatures with bare interlayer interactions, 
or at T = with gate screening (left). The value dAk/dfc|i at the zero crossing of as a function of temperature T with 
bare interlayer interactions or with gate screening (right). In both panels: d — as, Vs — 3, a = 0.3 and me/rrih ~ 0.07/0.30. 



superfluid states can be stable with momentum depen- 
dent interactioni22 We note in this context that in recent 
work, Sarma phases were found to be stable also in two- 
dimensional, two-band neutral Fermi systems i^i^ 



IV. SUMMARY 

We have studied a bilayer system of electron and hole 
layers spatially separated by an insulating barrier, where 
the electron and hole densities can be controlled indepen- 
dently and have analyzed s-wave pairing between elec- 
trons and holes as a function of average density and pop- 
ulation difference using mean-field theory. By solving 
the relevant energy gap equations we have compared the 
energy of the condensed phase called the Sarma phase 
with that of the normal state which is the sum of the 
electron and hole Fermi liquid energy described by the 
Hartree-Fock solution. We have included both inter- and 
intra-laycr interactions generalizing earlier work which 
did not include in-plane interactions.^'^ In this way the 
phase boundary for the ground state is established in the 
density-population polarization (rg — a) plane. The "lo- 
cal" stability of the Sarma phase was checked by calcu- 
lating the superfluid mass density, whereas the stability 
with respect to phase separation was assessed by calculat- 



ing the compressibility matrix. We have found that with 
bare Coulomb interactions the Sarma phase is always lo- 
cally stable due to a peculiar momentum structure of the 
gap function originating from the singular infrared be- 
havior of the Coulomb potential, and the simultaneous 
presence of a sharp Fermi surface at zero temperature. 

Employing a simple model of screening which intro- 
duces an infrared cut-off in the Coulomb interaction, we 
have found that some regions in the phase space be- 
come unstable. We interpret this as an instability to- 
wards an FFLO phase. Together with intra-layer in- 
teractions, the phase diagram in the crossover regime 
from the weakly interacting high density BCS limit to 
the strongly interacting EEC of dilute excitons has room 
for various phases. The topologically different SI and S2 
Sarma phases and FFLO are present with the inclusion 
of screening and intra-layer interactions. On the other 
hand, without any screening there is no instability to- 
wards FFLO and turning off intra-layer interactions the 
phase diagram does not show an S2 state. Currently, 
the experimental situation allows these systems to be 
realizedf^i*^!^ Quantitative comparison would require 
a more realistic model of screening, accounting for the 
condensed phase and finite width of the quantum wells, 
incorporating the disorder effects, and inclusion of spin 
degrees of freedom which may enter non-trivially when 
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FIG. 4: (color online) Phase diagram for inter-layer separation d = as with progressive refinement of the approximations. 
Superfluid (Sl/S2)-normal (N) phase boundaries are shown with red solid lines. A negative superfluid mass density showing 
a local instability is assumed to be towards an FFLO phase. SI, S2 and FFLO boundaries are shown with green dashed 
lines. The four panels corresponds to calculations including: bare inter-layer interactions only (upper left), bare intra- and 
inter-layer interaction (upper right), gate screened inter-layer interactions only (lower left) and gate screened intra- and inter- 
layer interaction (lower right). The a = line is special in the phase diagram and corresponds to the BCS state with equal 
populations, which is always stable. 



there are spin dependent interactions such as spin-orbit 
coupling. With the renewed mean- field phase diagram 
at hand, it would also be interesting to perform QMC 
simulations to probe the predicted phases. 
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